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Abstract. We give an approach for studying quantum gravity effects on black 
hole thermodynamics. This combines a quantum framework for gravitational 
collapse with quasi-local definitions of energy and surface gravity. Our arguments 
suggest that (i) the specific heat of a black hole becomes positive after a phase 
transition near the Planck scale, (ii) its entropy acquires a logarithmic correction, 
and (iii) the mass loss rate is modified such that Hawking radiation stops 
near the Planck scale. These results are due essentially to a realization of 
fundamental discreteness in quantum gravity, and are in this sense potentially 
theory independent. 

The first law of black hole mechanics is a result concerning a subspace of the 
solution space of Einstein's equations parametrized by a black hole's mass, charge, 
and angular momentum. Perturbations of these parameters are related, and this 
comprises the first law pQ. For the simplest Schwarzschild black hole the law is 

dM = — dA (1) 

Semiclassical arguments due to Bekenstein and Hawking [HE] have provided a physical 
and mathematical basis for identification of the surface gravity k with temperature 
and horizon area A with entropy. 

It is interesting to ask how the laws of black hole mechanics might emerge from 
a theory of quantum gravity. If such a theory is a conventional quantum system, one 
might expect that black holes are described by physical semiclassical states peaked 
at values of mass, charge, and angular momentum (M, Q, J). Such states may be 
denoted 

\M,Q,J) sc (t) (2) 

to indicate that they describe configurations peaked on classical values of the black 
hole's conserved charges, where t is a parameter characterizing the state's spread. It is 
reasonable to conjecture that the theory has operators corresponding to these classical 
observables, so that one can compute expectation values such as 

SC {M, Q, J\M\M, Q, J) sc = M + 0(M, A, t). (3) 
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Quantum corrections represented in the second term would in general be present 
because semiclassical states are not necessarily eigenstates of the relevant operators 
M, J, etc. These corrections would depend on t and on a fundamental discreteness 
scale A. Such states have been constructed for example for quantum cosmology 0]. 

What is required to provide a concrete implementation of these ideas is not just 
a quantum mechanical model of a static or stationary black hole, but a larger system 
capable of describing black hole formation and its subsequent evolution. One such 
system is the gravity-scalar field system in spherical symmetry, which is a non-trivial 
2— dimensional field theory. It has been well studied classically culminating in the 
discovery of scaling behavior at the onset of black hole formation [5] . One of the results 
is that black holes form with arbitrarily small mass for asymptotically flat initial data 
that is sufficiently massive. This is the so-called supercritical data. There is also a 
class of lower mass (subcritical) data for which there is no black hole formation. These 
are separated by a critical solution which is a naked singularity. It represents a finely 
tuned violation of the cosmic censorship hypothesis. 

Quantum gravity must modify these results due to at least two considerations: 
the existence of a fundamental length and a mechanism of singularity avoidance. 
Additionally these features should affect the thermodynamical properties of black holes 
in the Planck regime. This is based on the intuition that any mechanism for singularity 
avoidance must effectively be a form of gravitational repulsion. This effectively puts 
an upper bound on matter densities, similar to Fermi pressure [B]. This is turn should 
strongly affect the physics of the late stages of Hawking evaporation. 

Motivated by these considerations, we study here the effect of a fundamental 
discreteness scale on the thermodynamic properties of black holes. We show that the 
specific heat of a black hole can become positive as the temperature reaches the Planck 
regime, and that its entropy acquires logarithmic corrections. 

Our starting point is a Hamiltonian approach to the problem of quantum 
gravitational collapse for the same setting that has been well studied for classical 
collapse, namely Einstein gravity minimally coupled to a masslcss scalar field. A 
quantization of this system is available [7J, and this is what we use here. An important 
ingredient in the quantization is a class of operators that represent the radial function 

t) and its inverse in the spatial 3— metric 

ds 2 = q ab dx a dx b = dr 2 + R(r, tfdti 2 . (4) 

The full configuration space of the model is this radial field together with the scalar 
field <f>(r,t). 

Since the general case is dynamical, it is appropriate to consider thermodynamical 
variables associated with evolving horizons whose properties may be captured by 
gravitational phase space variables. It is possible to define quasi-local classical 
variables corresponding to energy and surface gravity that are known to satisfy "local 
laws" of black hole mechanics [S] . Using the quantization described in Ref . [7J it is 
possible to define operator analogs of these variables as a first step in formulating such 
laws in quantum gravity. 

A suitable definition of quasi-local energy for our purpose is the Misner-Sharp 
energy [3] 

E= ~R(l-q ab d a Rd b R). (5) 
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Other definitions, such as the Hawking energy, may also be used, but this does not 
effect the approach we take in this paper. This expression together with the definition 

k = l -q ab d a d b R (6) 

for the surface gravity of a trapping horizon leads to the quasi-local first law [8] 

dE=^-dA- l -pdV. (7) 

Each quantity here is evaluated on the trapping horizon; A is its area, V = 4/3irR 3 
and p = \j a \ the norm of the matter current (j a = q ab db4> for the scalar field). 

For the case where the trapping horizon is slowly varying or not at all, the second 
term in Q vanishes, and k becomes inversely proportional to the horizon radius Rh- 
This is familiar from the Schwarzschild solution. Although our approach applies to 
the general dynamical case, for simplicity we consider here the situation where the 
radial size of the local horizon is approximately constant. Thus on the horizon we 
have 

E= l -R H , K= R 1 h - (8) 

With these consideration in mind, we propose the following state dependent 
definitions of trapping horizon entropy S and specific heat C entirely within the 
quantum theory. 

ds = Sl, c = §§, (9) 

(k) d{K) 

where E and k are suitably defined operators representing the corresponding classical 
quantities. The differentials in these equations are to be taken with respect to the 
horizon size parameter characterizing the state. This could be the parameter at which 
a semiclassical state is peaked, or an eigenvalue in case of a state in which the operators 
are diagonal. We shall use the latter to show that the familiar entropy formula modified 
by logarithmic corrections is recovered, and that the specific heat exhibits a phase 
transition from negative to positive values near the Planck scale. 

In the quantization described in [7] , basis states are associated with a radial lattice 
of points {rfc} (k = 1 • • • N), with field values at the points prescribed. For a theory 
with only one field degree of freedom, the basis states are 

\R 1 ,R 2 ,---Rn)- (10) 

For matter coupling, the basis is simply a product of states of this type. 

The operators representing configuration variables such as R are diagonal on these 
states 

R(rk)\Ri, R2, ■ ■ ■ Rn) — Rk \Ri-,R2, ■ ■ ■ Rn), (11) 
and inverse configuration operators are realized as 



\Ri, i? 2 , • • • Rn) 



{\R k + X\^ l -\R k -X\^ 1 



n «/(*-!) 

\Ri,R 2 , - ■ ■ Rn) 



I 

2A 

/(i2jfe,A,0 ]i2i, J Ra,---J2 JV ) (12) 
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Figure 1. The specific heat as a function of horizon radius R/X for the case 
I = 1/5. The singularity and positive region for small horizons is apparent. For 
larger horizons the function matches the negative semiclassical result. 



There are two parameters in the definition of the 1/ R operator. One is the fundamental 
discreteness scale A and the other is a number < I < 1 that characterises an 
ambiguity in defining this operator. For r >> A, the function / behaves like but 
for sufficiently small R it has I and A dependence. This is illustrated by the expansion 
for R » A : 

^• A -'>4 +i? w^ +o < A4 » (i3 > 

The entropy and specific heat for a trapping horizon of radius Rh can be 
computed using equations ([9]). The quantum gravity modifications in these quantities 
may be illustrated by considering a lattice basis state \Ri, R2, ■ ■ ■ Rn), where the field 
value Rr occurs at one of the lattice points. We find 

,0 = d(E) = c 2 R H dR H 

(X/R) 2G l + (2J-l)A 2 /6fi&' 1 1 

— 2G H ~ UP 1uRh - (15) 

This shows that the entropy has a In Ah correction that depends on the discreteness 
parameter appearing in the inverse R operator, as well as on the / ambiguity in the 
definition of this operator. This correction vanishes for I = 1/2. Although logarithmic 
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corrections to horizon entropy are known in various approaches [10] . ours differs in 
that its origin is intimately tied to the l/R operator, which in turn is connected to 
singularity avoidance [7J. 

The proposal ([9| for the specific heat C leads to the quantum gravity modified 
formula 

r = 2r4h1 na\ 

2R 2 H l + {21 - l)X> ( ' 



This formula ( 16 1 is derived from an expansion of the function f(R, A, I) in powers of 
A rather than the full expression. 

For large black holes Rh » A this gives the familiar negative specific heat, 
independent of A and I. For I > 1/2, C is finite and negative for any horizon size. The 
most dramatic case (illustrated in the figure) occurs for I < 1/2. The specific heat 
becomes singular when a black hole reaches the size 



1 — 21 

Rait = A V^~' (17) 

and becomes positive as its size reduces further still. This suggests a second order 
phase transition in quantum gravity to a phase where a black hole becomes a 
thermodynamically stable object. 

If the exact formula for f(R, A, I) is used rather than its expansion in powers of A, 
we find that the specific heat has the same qualitative behaviour as in Fig. [I] but with 
the two important differences: the specific heat is negative for large R and positive 
for small R for all values of the parameter I (not just for I < 1/2), and the location of 
the transition is closer to R = 0. 

A further consequence of these arguments is a modification of the black hole mass 
loss rate. The usual formula may be written as 

dt t P \ M J y ' 

where a is a positive dimensionless constant and Mp,t p are the Planck mass and 



time. The eigenvalue of the l/R operator (12) suggests the following modification of 
the mass loss formula: 

-i 2 



dM _ M P 
dt tp 



1 (1 - 21) 



Mp \m--W2W)\ ■ (19) 
It is apparent in this modification that the mass loss rate vanishes for sufficiently small 
R for I < 1/2, and that the usual formula is recovered for large R. We note again 
that if the exact formula is used the rate vanishes at M — rather than at a value of 
the order of Mp. In either case the evaporation rate slows to zero and a "black hole 
explosion" is averted. 

Although our arguments are exhibited using basis states in a polymer quantization 
of the spherically symmetric sector, it is expected from considerations of the inverse 
scale factor in quantum cosmology that qualitatively similar results will hold for 
semiclassical states jTTJ [T^] . Thus these results may reflect a true feature of quantum 
gravity. 

Our considerations are independent of the value of the Planck scale and therefore 
apply also to large extra dimension settings. Indeed they depend fundamentally on 
only two factors: the existence of a fundamental discreteness scale and its effect on 
the realization of operators such as l/R. It is in this sense that our arguments are 
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potentially theory independent. The phase transition evident in the specific heat has 
potentially important consequences for black hole production at the Large Hadron 
Collider - it suggests that Hawking evaporation will cease after black holes reduce 
to the critical size given by the above formula. Furthermore it provides a quantum 
gravity basis for the possibility that missing energies at LHC collisions may be carried 
by Planck scale remnants [13j . 

It is of interest to contrast our results and formalism to other works that take 
a different approach to the problem of black holes in quantum gravity. Some early 
works considered computing black hole temperature corrections due to back reaction 
of Hawking radiation on the Schwarzschild geometry. This is then used to compute the 
specific heat, which is found to exhibit a divergence and sign change at a critical value 
of temperature [14] . This result is qualitatively similar to what we report, although 
our approach involves non-perturbative quantum gravity considerations, rather than 
a perturbation theory calculation of back reaction. 

Other ideas have involved the use of modified uncertainty relations to argue 
for an upper bound on black hole temperature, in some cases motivated by non- 
commutative geometry 15J. In these approaches a model of non-commutative gravity 
is used to extract an analog of a Schwarzschild solution. This solution is then used 
to extract thermodynamical properties of black holes is the usual manner. Again 
some of the qualitative features of specific heat appear similar to what we obtain here 
from canonical quantization of the spherical symmetric sector. These similarities of 
our results with other significantly different approaches suggest that some universal 
considerations involving a fundamental length scale are at play. This is somewhat 
reminiscent of the entropy-area relation for black holes in that whatever postulates 
are used for microstates, the entropy turns out to be proportional to area|16j. 
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